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We demonstrate the realization of nearly massless electrons in the most widely used device material, silicon, 
at the interface with a metal film. Using angle-resolved photoemission, we found that the surface band of a 
monolayer lead film drives a hole band of the Si inversion layer formed at the interface with the film to have 
nearly linear dispersion with an effective mass about 20 times lighter than bulk Si and comparable to graphene. 
The reduction of mass can be accounted for by repulsive interaction between neighboring bands of the metal 
film and Si substrate. Our result suggests a promising way to take advantage of massless carriers in silicon-based 
thin-film devices, which can also be applied for various other semiconductor devices. 
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The ultimate performance of electronic devices is largely 
governed by the effective mass of charge caiTiers. There- 
fore, the never-ending quest for higher performance devices 
has looked after high speed carriers with light effective mass 
Uj] . This quest may have reached its ultimate goal through the 
recent findings of massless electrons with their speed close to 
light in graphene t^lM and bismuth compounds L4-7J. How- 
ever, for these materials various difficulties exist in promptly 
realizing practical and mass-producible devices. 

In general, the effective mass of electrons is determined by 
their energy dispersion. Especially, the edge structures of con- 
duction and valence bands near the energy gap (Eg) are im- 
portant. In the band-edge region, the dispersion is sensitive to 
the interaction between the bottom conduction and the top va- 
lence band as described well by the k ■ p theory |l8|]. This stan- 
dard theory provides an approximate analytic expression for 
band dispersion under the interband interaction around high- 
symmetry points of momentum (fc) space |8]. The interaction 
energy, ±^/{Eg/2)^ + (P • k)^ (P, the optical matrix ele- 
ment), gives rise to repulsion between bands and provides the 
linear component in dispersion [8]. As Eg approaches zero, 
the dispersion becomes nearly linear E(k) Ri P ■ k with a much 
lighter effective mass of electrons m* = h^Eg/{2P'^). This ef- 
fect explains the trend that a semiconductor with a narrower 
gap has a lighter effective mass |9] and the zero-gapsemicon- 
ductor of graphene has a negligible effective mass ^\M- This 
theory also tells us that if the energy gap could be controlled, 
the effective mass would be tuned through the interband inter- 
action. While the energy gap of bulk materials cannot be tuned 
easily as fixed by their crystalline structure, here we show that 
a similar effect can be obtained at a semiconductor interface 
where a proper interface state is formed within the band gap 

iMl. 

We chose an ultrathin (only single layer) Pb film on n-type 
Si(lll) substrate, which has ideally abrupt interfacial struc- 
ture and strongly dispersing metallic electron bands within 
the Si band gap lll3l4l5ll . For this system, the film-substrate 
interaction, which is exploited in the present work, was al- 
ready invoked to explain the anomalous superconductivity at 
the two-dimensional (2D) limit jla llTll . Figure [TJ a) illus- 
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FIG. 1: (color online), (a) Crystal- and (b) reciprocal- lattice struc- 
tures for Pb/Si(l 1 1) at 1.2 ML with a ^ x ^ symmetry 1 13]. Grey 
(dashed) line represents the \/7x -\/3 (1 x 1) surface unit cell. ARPES 
data collected along (c) long and (d) short arrows in (b) crossing 1x1 
(f) and y/l X \/3 (Fs) Brillouin zone centers, respectively. The 
dashed line in (c) denotes the Si valence-band edge. A-shaped band 
with a linear dispersion on each Fg in (d) is labeled R and two other 
gap-state bands in between are denoted by SI and S2. Both ARPES 
images are symmetrized with respect to the origin, (e) Band bend- 
ing in the present n-type Si substrate calculated by the measured Ef 
position and the bulk doping concentration (231 . 



trates the atomic structure of the interface with the Pb den- 
sity of 1.2 monolayer (ML) (7.84 atoms/nm^). Pb atoms are 
densely packed within a single layer upon the bulk-terminated 
Si(lll) surface 111311 . Most of Pb atoms sit on top of underly- 
ing Si atoms (Tl or Tl' sites) while part of them are slightly 
displaced (Tl') by additional Pb atoms in hollow sites (H3), 
which lead to the formation of a uniform ^/l x Vs unit cell 
(grey lines) f 18!]. Due to the anisotropy of this unit cell and the 
three-fold symmetry of the substrate, one inevitably obtains 
triply rotated domains il5Lll8ll . which should be considered 
in interpreting experimental data. 
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FIG. 2: (color online), (a) ARPES data of the R band. Data below 
the dashed line are those symmetrized with respect to Ts to elimi- 
nate the strong neighboring feature S2. (b) Spectral peak positions of 
the R band (open circles) extracted from MDCs shown in (c). Data 
points given in closed circles are those symmetrized. The dispersion 
of a normal Si LH band (black parabola) is compared to show the re- 
duction of curvature. Red and blue lines are the results of k-p model 
fit I2al . The error bar in k is less than ±0.02 A"'^, which brackets the 
quantified effective mass (orange lines) (see also the supplementary 
information L18.1 ). Red lines overlaid in (c) show fits by Lorenzian 
functions. 



Angle-resolved photoemission (ARPES) measurements 
were conducted in an ultra-high vacuum chamber (6.5 x 10^^ 
torr) equipped with a hemispherical electron analyzer (R4000, 
VG Scienta) and a high-flux He discharge lamp for 21.2-eV 
photons. The samples were cryogenically cooled down to 90- 
100 K for measurements. The overall energy and momentum 
resolutions were better than 20 meV and 0.02 A~^. Figure 
[He) shows band dispersions along the high-symmetry direc- 
tion through the center of Brillouin zone f [the long arrow in 
Fig. [lib)] measured by ARPES. Near f , an intense parabolic 
band is readily found, which is the well-known direct tran- 
sition from bulk Si valence bands il9ll . Another feature with 
strong intensity appears near each zone boundary (A') with lit- 
tle dispersion around 0.7 eV, which is the covalent-type bond- 
ing state between Pb and Si lll4ll . In contrast, around K, 
there are parabolic bands dispersing toward the Fermi energy 
(Ep), which were identified as due to 2D metallic electrons 



localized within the Pb layer (called S2 hereafter) Ill5ll . This 
metallic surface state induces a huge upward band bending in 
the n-type substrate, and forms a p-type inversion layer (IL) 
as shown in Fig. [TJe) 12011 . The strong band bending yields 
a triangular potential well to confine electrons in the inter- 
face IL into 2D quantum well states (QWS) EESl- An IL 
state at interfacial layers is typically not easily probed by very 
surface sensitive ARPES but, for an ultrathin metal film like 
single-layer In on Si(l 1 1), a previous ARPES study traced its 
detailed dispersion |21]. We also found such IL states; two 
hole-like bands are observed with weak intensity between the 
bulk and surface bands in Fig. [TJc). A relatively prominent 
band is identified as the so called heavy-hole (HH) bands 12111 
but the other branch with a lighter effective mass, light-hole 
(LH) band, is barely observable as a weak and broad feature 
on top of the Si bulk band. 

Although the LH band is not clear enough here [Fig. [TJc)] 
due to its weak intensity, we can map its dispersion clearly 



at the center of the surface Brillioun zorie (Ts) away from T 
as translated by the surface periodicity il2l 12111 . Figure [TJd) 
shows band dispersions taken through two such Fg's and the 
M point [the short arrow in Fig. [lib)]. There are strongly dis- 
persing bands with dominant intensity, crossing Ep at ±0.21 
A^^, which correspond to the S2 state for Pb metallic elec- 
trons. In addition, two other bands are identified, SI folded 
back near Ep with respect to M and R folded with respect to 
Ts- Surprisingly, their dispersions are A-shaped and appar- 
ently very linear, especially for R. The detailed dispersion of 
R, located on Ts and related to the Si LH band below, can 
be shown more clearly and quantified by the peak positions of 
momentum distribution curves (MDCs) in Fig. I2c). The dis- 
persion is indeed linear within the experimental uncertainty 
as shown in Fig. Wih) (open circles). Thus, the effective mass 
value would be extremely small and cannot be quantified by a 
simple parabolic fit of the dispersion as in the case of graphene 
yfl and Bi compounds ||5|-(7|]. Instead, the Fermi velocity can 
straightforwardly be extracted from the slope of the band and 
is as high as 4.6 ± 0.4 x 10^ m/s, which reaches to half of 
those in grapheneO, [l8[] and BinqSbo.i 15|], and is similar to 
those in Bi2Se3 |6] and Bi2Te3 [7]. The effective mass value 
of this linear band will be discussed further below. 

Since the R band with a striking dispersion cannot be found 
in Si bulk crystals [Fig. Ob)] and on clean Si(lll) surfaces 
l,19il . it must be due to electron states localized at the interface 
(Si ILs) or within the Pb layer To identify the origin of each 
band further, we performed theoretical calculations based on 
density functional theory (DFT). The calculations were per- 
formed using the Vienna ab-initio simulation package within 
the generalized gradient approximation and the ultrasoft pseu- 
dopotential scheme with a plane wave basis [22J. The surface 
is modeled by a periodic slab geometry with 6 and 12 Si layers 
whose bottoms are terminated by hydrogen. 

Figure [3 a) shows the calculated constant energy contours 
at Ep. There are two kinds of surface-state Fermi contours, 
point-like crossings and wavy lines repeated following the sur- 
face periodicity. Since ARPES detects signals coming from 
all triply-rotated domains, the expected Fermi contours result 
in a complex pattern [Fig. [3jb)]. Nevertheless, the agree- 
ment between the calculated and experimental Fermi contours 
is excellent in Figs. [3jc) and [3d). The point-like crossings 
and the details of the wavy contours are all precisely repro- 
duced, which correspond to S 1 and S2, respectively. The dis- 
persions of S 1 and S2 are also reproduced quantitatively well 
in the calculation [Fig. [3f)]. These states come from the Pb 
layer due to in-plane Pb 5p orbitals. In fact, SI and S2 bands 
are doubly-degenerated with contributions from domains B 
and C overlapped exactly [Fig. [3jf) and green and blue ones 
in Fig. [3jb)]. The contribution from domain A, which pass 
through two f 5 points [Fig. [3e) and red ones in Fig. [3b)], is 
not degenerated. For this domain, the Pb-derived band (S3) is 
located just above Ep and two Si-derived bands below Ep at 
Ts [Fig. [3e)]. The calculated electron density for the latter 
is distributed within Si layers and their dispersions are consis- 
tent with the Si IL states (LH and HH) mentioned above. In 
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FIG. 3: (color online). Constant energy contours at Ef calculated 
from DFT for the (a) single-domain and (b) triply-domain surfaces. 
Contours from each domain (A, B, and C) are indicated by different 
colors in (b). (c) ARPES constant energy contours at Ef- The DFT 
contours from (b) are superimposed. The raw data (region without 
DFT contours) are symmetrized reflecting the fundamental mirror 
symmetry, (d) Enlarged data around M for the detailed comparison 
between theory (while lines) and experiment (red circles extracted by 
MDC peak positions). DFT band dispersions along the arrow in (b) 
for (e) domain A and (f) domains B and C, which are overlapped 
together in (g). The shaded area represents the projected bulk band. 
For the best match with the experimental data the calculated surface 
state energies are rigidly shifted by 120 meV with respect to Si bands. 



particular, one with a sharper dispersion (LH) is very similar 
to R as compared in Fig. |2b). The calculation, thus, clearly 
identifies SI, S2, S3 and R bands as three 5p states of the Pb 
layer and a Si-derived state, and reproduces the dispersions of 
SI, S2, and S3 quantitatively well. This result is consistent 
with the very recent work in Ref. 11231] for the single-domain 
surface. 

One notable feature of the calculation is that the bottom of 
the S3 band lies on Ep, namely, just above the experimentally 
measured band R with a very small energy separation. Thus, 
S3 has a chance to interact strongly with R from Si subsurface 
layers through the interband interaction |8]. As mentioned 
above, this interband interaction determines the effective mass 
of the bands according to the k ■ p model [8]. Indeed, such an 
interaction between the band of a Ag monolayer film and the 
band of a Si(l 1 1) substrate was observed recently and affects 
the dispersion of a S3-like surface state of Ag ll24ll . 



The fitting of our experimental dispersion using the k ■ p 
model 1 25] with the band-gap size determined by the calcu- 
lated energy position of S3, 23 meV, reproduces well the lin- 
ear dispersion of the R band as shown in Fig. |2|b), yielding an 
extremely light effective mass of 0.0074 ± O.OOlSmg, where 
TTie is the electron rest mass. In fact, the k ■ p fitting is simple 
enough to treat the band-gap size as another free parameter 
and the band gap is reliably determined to be within 18-24 
meV il8|] in good accord to the theoretical value. The effec- 
tive mass determined is 20 times lighter than the normal Si LH 
band (O.lSmg) 19]. Note also that, for the HH bands, the in- 
terband interaction is forbidden from the band symmetry i^. 
Therefore, we conclude that the R band with an extremely lin- 
ear dispersion is the modified LH band through the interband 
interaction with S3, while we do not know why the HH band 
itself is not observed around Ts- 

The k ■ p method has been widely applied to various 
semiconductors and the smallest effective mass reported is 
0.016me for InSb Il9l[l8|]. This yields the carrier mobility 
as high as 30,000 cm^ V~^ s~^, about 2 orders of magnitude 
higher than bulk Si uJ]. The present effective mass for the Si 
IL, thus, indicates the possibility of an ultrafast carriers at Si 
interfaces. This method has not been applied to graphene and 
Bi compounds, for which only the effective mass values mea- 
sured by transport experiments are available; about (0.002- 
0.009)me in a similar range with the present value ],2. ,5> , 1 8|] . 
In contrast, the DFT calculation shown above does not repro- 
duce exactly the linear dispersion of the LH band, yielding an 
effective mass of roughly O.Qlme- We think this discrepancy 
is due to the limitation of the present slab calculation with 
only few Si layers in reproducing properly the QWS of the 
much thicker IL. 

The proposed mechanism for the linear dispersion is sum- 
marized in Fig. SJa). The contact of a metal film to a semi- 
conductor substrate plays two key roles; (i) to induce a strong 
IL near the interface lifting the hole bands of the substrate 
up close to Ep and (ii) to provide a proper band (S3 here) 
within the band gap inducing the strong repulsive interaction 
with hole bands for linear dispersion. In order to verify this 
mechanism, we took advantage of diverse surface band struc- 
tures in metal/Si(lll) systems. We examined three repre- 
sentative metal overlayers — Pb, Au, and Ag on Si(lll) — 
with a common \/3 x \/3 symmetry but with different sur- 
face bands. Figures |4|b)-0td) describe corresponding band 
dispersions around F taken from previous reports in literature 
llHllllIS]^ The a/3 X V^-Pb system (with the Pb density 
of 4/3 IVIL) has no bands close to Si LH bands ],15.] while the 
others (at 1 ML) have adequate bands from the films similar 
to the present case to allow the interaction with the LH band 
(band gaps of 40-60 meV) 12611271] . The consequence of such 
a difference on the dispersion of the LH band is remarkable. 
In Fig. life), our ARPES data for the Pb system show a nor- 
mal parabolic dispersion consistent with bulk Si [8]. However, 
for Au and Ag systems, dispersions are obviously much more 
linear [Figs. Uf) andUg)]. These results not only provide 
solid evidence on the validity of the proposed mechanism but 
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FIG. 4: (color online), (a) Schematic illustration of the mecha- 
nism for a linear hole band dispersion in semiconductor ILs. (b)-(d) 
Schematics of surface (SS) and LH band dispersions near F in Pb, 
Au, and As single-layer films on Si(l 1 1) with a common y/S x \/3 
symmetry lla. 1241 1291 . (e)-(g) Corresponding ARPES dispersions 
of LH bands measured along the long arrow in Fig. [Hb) across T. 



also suggest the possibility of applications for a wider range 
of materials. 

Our result suggests an unprecedented way to ultrafast elec- 
tronic devices based on Si. Since this technique directly con- 
trols the carrier mass at the interface by a metal overlayer, it 
does not require any artificial modification or chemical engi- 
neering on the Si bulk lattice itself. Furthermore, this mecha- 
nism, in principle, can be generally applied to various other 
semiconductors even with thicker films llllll . provided that 
proper interface states exist. For example, this may solve the 
notorious imbalance between the hole and electron mobility 
of III- V compounds [1]. 
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